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Abstract 

We study numerically the induction mechanisms generated from an array of helical motions 
distributed along a cylinder. Our flow is a very idealized geometry of the columnar structure that 
has been proposed for the convective motion inside the Earth core. Using an analytically prescribed 
flow, we apply a recently introduced iterative numerical scheme [TB] to solve the induction equation 
and analyze the flow response to externally applied fields with simple geometries (azimuthal, radial, 
etc.). Symmetry properties allow us to build selected induction modes whose interactions lead to 
dynamo mechanisms. Using an induction operator formalism, we show how dipole and quadrupole 
dynamos can be envisioned from such motions. 
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I. INTRODUCTION 



The understanding of the self-sustained dynamo of the Earth is still a major challenge. 
Following Larmor's original hypothesis pQ, it is supposed to originate in the convective mo- 
tion inside the liquid iron core. There, the induction due to fluid motions may overcome the 
Joule dissipation and a dynamo can be generated. Although knowledge of the fluid motion 
is essential because it drives the magnetic induction, the structure of the core convective 
motion is not precisely known. It is due to the extreme range of parameters in this problem: 
the Earth is in rapid rotation, and the electrical conductivity of molten iron is large, but 
finite. As a result, two important dimensionless parameters of the problem are very small: 
the Eckman number may be as low as 10~^^ and magnetic Prandtl number of the order 
of 10~^ [2] are usually quoted. In addition, the source of the convective motion is mixed 
and still debated: thermal and compositional convection contribute, with heat exchanged 
at the boundaries (solidification of the iron at the solid inner core, transfer at the mantle 
boundary) and heat is also released in the bulk because of radioactivity [3]. The Rayleigh 
and Nusselt numbers are not precisely known, and the coupling with the magnetic field can 
result in significant changes in the structure of the convective flow. However, many models 
of the Earth dynamo have invoked the columnar flow structure first derived by F. Busse |1] 
at onset of convection in a spherical Couette geometry These columns may persist for 
a system moderately above the onset of convection ^ [7] although their existence in the 
actual Earth core is still an open question [S]. There are also arguments for a subcritical 
dynamo bifurcation [9J, in which the columnar structure may be stabilized by the large scale 
magnetic field [TO] . 

The focus of this paper is on the induction mechanisms that result from the existence 
of such a columnar structure. We consider a highly idealized situation in which columns 
with helical motion are distributed along a cylinder, as helicity is known to favor dynamo 
action [TT] . The geometry is chosen to be cylindrical so that the conical and curves conditions 
of the Earth are ignored. We also do not take into account the drift of Busse's rolls, nor the 
presence of a zonal wind. The flow is prescribed, and we analyze its response to externally 
applied fields with simple geometries (azimuthal, radial, etc.) in order to understand the 
induction mechanisms in this system. For instance, several numerical simulations (e.g. [12] 
and references therein) have noticed that both dipole and quadrupole dynamos are possible, 
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with thresholds in close proximity. Note that the periodic array of helical columns in the 
Roberts flow [13] and Karlsruhe dynamo [Hj cannot generate an axial dipole. We show here 
that the annular geometry can lead to both axial dipoles and quadrupoles. We describe in 
detail the flow structure and the numerical method in section II. In sections III we present 
our results concerning induction responses to simple applied flelds. Section IV is devoted to 
the study of possible dynamo action in this type of flows, based on an operator formalization. 
Conclusions and possible extension of the study are presented in section V. 



II. SYSTEM AND METHODS 



A. System 

The flow geometry is made of a system of cylindrical columns forming an annulus, as 
shown in flgure 1. The columns are grouped in pairs of cyclonic and anti-cyclonic roll, for 
which the axial flow is reversed. The velocity is assumed to be stationary and is expressed as 
the sum of a contribution due to the circular motion of the fluid in a column (the rotational 
component V^(r)) and of a contribution due to an axza/ motion, which we label V'^(r). In 
an Earth-like geometry, V'^(r) would be generated by the Eckman pumping at upper and 
lower boundaries. One thus writes 



V(r) = V V) + V) , (1) 

We chose analytical expressions of the flelds such that the components are separately 
divergence-free, V ■ V^(r) = V ■ V"^(r) = 0. The coefficient ^ measures the intensity of the 
axial motion compared to the rotational one. 

The rotational part in columns of height 2H is expressed in cylindrical polar coordinates 

as: 
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where n is the number of column pairs, R is the outer radius of the annulus, d is the thickness 
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of the region in which the columns are confined (within radial distances between R — d and 
R). The velocity is set to zero outside the domain R — d < r < R. 

For the axial flow, we consider two cases. In the first one, the columns have a height 
2H and the axial flow has a defined direction within each column and reverses between 
neighboring columns. It corresponds to the geometry sketched in figure [l|^a). Figure [l](c) 
and (d) show respectively a cut of the flow in the planes z = and ^ = 0. The helicity 
(7i = V ■ (V X V)) has the same sign in each column, negative in this case (a column with 
positive axial velocity rotates with negative vorticity). The corresponding poloidal velocity 
is given by: 



In the second case (sketched in figure [T](b)), the axial flow is reversed about the plane z = 0, 
hence having symmetries similar to Busse's columns in rotating convection. Figure [l](e) 
shows a cut of the flow in the plane 6 = 0. The helicity is therefore negative in the columns 
of the upper half of the cylinder and positive in the lower half. The axial velocity is then 
written: 



first choice of the axial velocity field, and Tl the configuration obtained with the second 
choice. We shall look for stationary solutions of the induction equation: 



where rj is the magnetic diffusivity of the fluid. The boundary conditions are such that 
the medium inside the annulus has the same electrical conductivity as the fluid while the 
outside medium is insulating. The magnetic permeability is equal to that of vacuum in the 
whole space. 




(3) 



<^ V.^'iz, 6) = -—sm{n9) cos j (4) 
r/(z,^) = cos(n^) sin(^^) 

In our study, we call T| the configuration obtained with the rotational velocity and the 




= = V X (V X B) + r/AB , 



(5) 
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We stress that once the radius R of the cyhnder and the aspect ratio H/R = 2 are 
fixed there are still many independent parameters which may be varied: the number 2n of 
columns, their aspect ratio d/R, the magnitude of the axial flow compared to the rotational 
one, etc. We shall concentrate on the geometry portrayed in flgure 1: four pairs of columns 
with relative thickness equal to 0.4 (yielding a square aspect ratio for the columns), and 
^ = 1.25 (this sets the ratio of the maxima of axial to rotational velocities to 0.7, close 
to the values in the existing experimental dynamos [HI [151 123]). The only remaining free 
parameter is the amplitude of the velocity fleld, which is non-dimensionalized in the form 
of the magnetic Reynolds number = V^.^^R/ri. l^ax is the maximum velocity in the 
domain of the flow, and we write V = VmaxV. 

B. Iterative procedure 

We use here the iterative technique introduced in [T7] . The reader is referred to it 
for a detailed presentation, and for an evaluation of its performance compared to standard 
analysis in magnetohydrodynamics, such as the linear stability analysis of the induction 
equation. We stress that the technique is best suited for the study of stationary solutions, 
and only recall here its basic principles. 

We consider the response of the flow to an applied magnetic fleld Bq, looking for the 
induced fleld B which solves the stationary induction equation 

V X (V X B) + 77AB = -V X (V X Bo) . (6) 

The result is expressed as the integer series 

00 

B = J]B, with |Bfc| ~ 0(R^)So , (7) 

k=l 

The contributions B^ are computed iteratively from a hierarchy of non-dimensional Poisson 
equations 

ABfc+i = -R^V X (v X Bfc) , A; > (8) 
which can be solved for any given set of boundary conditions. 
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C. Numerical simulation 



Applying a standard Poisson solver to equation would require to write the complete 
set of conditions for the magnetic field at the boundaries of the computational domain. In 
practice, the boundary conditions are more readily expressed in terms of electric currents 
and potentials. We thus implement the following sequence: 

(i) The electromotive force (e.m.f. in units of Knax-Bo) induced by the flow motion is com- 
puted from e^+i = v x B^. 

(ii) Electric current being divergence free, the distribution of electric potential is obtained 
from A0fc_|_i = V ■ (v X Bj^), with von Neuman boundary conditions (n ■ V0 = n ■ (v x B), 
with n the outgoing normal of the domain). 

(iii) Induced currents (non-dimensionalized) are then computed as from Ohm's law 
jfc+i = — V0A:+i + Gfc+i, and used to compute the magnetic field B^+i from Biot and Savart 
law. 

All calculations are made with a magnetic Reynolds number equal to unity (we set i? = 1, 
Vmax = 1 and ?7 = 1). Its actual value enters only in the final step, when contributions are 
collected and the integer series is computed: B(Rm) = J2k^k(Rm = l)Rm- This approach 
requires that the series converges, and sets an upper value for the magnetic Reynolds num- 
ber Rjn (Rm < R-m)- Wc havc fouud that the radius of convergence Rm* is of the order 
of 30; for higher magnetic Reynolds number values we have shown in that Fade ap- 
proximants [18j still give results in remarkable agreement with the solution of the induction 
equation computed without approximation. 

III. INDUCTION PROCESSES. 

A. Induction from a toroidal applied field. 

When an external magnetic field is applied in the azimuthal direction (Bo = -Bo^e), one 
expects the generation of an azimuthal current, in very much the same manner as the a 
effect [in] operates in the Roberts flow [13] and in the Karlsruhe dynamo [HI ^7\. This is 
because there is no conceptual difference between a horizontal field transverse to the columns 
in the Karlsruhe geometry and a toroidal field applied in the geometry considered here. In 
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the azimuthal direction, there is a scale separation between the size of one column and the 
circumference of the annulus, so the results of mean-field MHD, as for instance explored 
in [20], should apply. 



1. Induction at first order 

We start with induction at first-order, i.e. the solution of the first iterative step of the 
induction equation: 

ABi = -V X (v X Bo) (9) 

The topology of Bq is chosen identical to the one that would be generated by the flow of 
an electric current in a rod of radius a located along the z-axis 

{ -ee {r <a = 0.1) 
Bo = i?o g ^ ^ (10) 

-ee (a < r) 
V r 

The induced field Bi is shown in figure [2|a) and (b). One observes two major effects: in 
the axial direction the field lines are stretched by the vertical velocity gradients, and in the 
transverse plane they are advected by the rotational motion in the columns. The mechanisms 
are schematically shown in figure [2](c) and (d); Bi^z (figure [2|^c)) results from the action of 
the axial flow, as can be seen by writing the axial projection of equation (|9]): 

ABi,, = -Bo%vf . (11) 

Note that since the stretching is largest at the boundary between two columns, the induced 
field is also maximum at these points, therefore out of phase by a quarter of a period with 
respect to the velocity field. 

For the transverse components of the induced field, (r, 6) directions, all azimuthal gradi- 
ents of the velocity contribute a priori 

(ABi)(,,e) = Bo^de^rf.^e) + ■ VBq - Bo^de^rf^^e) + " deB, . (12) 

However, the picture is simpler when restricted to induction in the z = plane. In this case 
the last two terms of the above equation do not contribute because v"^ is an odd function 
with respect to z. In particular, one obtains for the field induced in the radial direction 

= {dev^ - vf) + ^5o^ = B^^dev^ , (13) 
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which shows that (i) the stretching by the azimuthal gradient of is compensated by the 
advection of the apphed field; (ii) the net effect is due to the azimuthal modulation of the 
radial flow (figure |2](d)). At first order, the rotational flow component generates in the 
center of the columns an induced field perpendicular to the applied one, a mechanism which 
is quite helpful in the understanding of field expulsion by vortices. 



2. Second order effects 

We now turn to the analysis of the induction at second order of iteration, i.e. the 
induction resulting from the presence of Bi. The electromotive force e2 = v x Bi is shown 
in figure [3]^a,b). At this stage, the induction from each column add- up cooperatively. An 
azimuthal current J2,e is induced - figure [3|^b) - parallel to the applied field Bq. It has 
the symmetry of the flow helicity. As we have verified, it is reversed as either v"^ or 
is reversed but not when both change sign. J2fi generates the axial field -82,2 shown in 
figure [Sj^c). This induction by helical motion, schematized in figure |3|^e), was first proposed 
by Parker [21], and evidenced in the VKS experiment [22]. It is also in agreement with the 
a effect introduced in the framework of mean-field magnetohydrodynamics [19], for the 2D 
array of columns in the Roberts flow [13J, and for the columnar ring discussed here [20]- 
Following this analogy, we will call this effect a-effect in the rest of our study. 

The electromotive force e2 also generates axial currents which create the azimuthal field 
B2fl shown in figure [sj^d): it is largest in the center of the columns and its direction is opposed 
to the applied field. Its generation - schematized in [3|^f) - shows that i?2,e is opposed to 
the applied field even when the rotation of the columns is reversed, or when the axial flow is 
reversed or even suppressed. The effect traces back to the expulsion of magnetic field lines 
by vortical motions, as shown by numerical [21] and experimental ^25j studies on an isolated 
vortex. It is a second order effect: as pointed out previously, the rotational motion of the 
column can be thought of as turning the applied field lines by 90 degrees at each step, so 
that after two iterations the induced field is directly opposite to the applied one. 

We stress that the response of the flow to an applied field is more complex than the 
generation of an azimuthal current by an 'a' effect due to the helical motion in the 
columns; this effect is indeed present, but the net current J2 is related to the applied field 
by a full second rank tensor, 3 2 = c"[a]Bo. In the tensor [a] not all components are due to 
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the flow helicity and some components correspond to an expulsion effect. 



3. Azimuthal averages 

In the remaining of our study, we consider the contributions to the induction averaged 
over the azimuthal angle 6, as they present a somewhat clearer picture. Figure |4] shows the 
averaged order 2 fleld, deflned as: 

{B,)ir,z)^^ j\eB,{r,e,z) (14) 

One clearly sees the induced azimuthal current {J2,e) and the associated magnetic fleld with 
a dipolar structure. In a flnite geometry with insulating surroundings, the currents (J2) 
must remain within the flow volume, leading to poloidal current loops (arrows in ^a)) 
which participate to the expulsion of the applied fleld. 

In order to quantify the importance of each effect, we have plotted in figure [5]the evolution 
of the components of (B2) with the radial and axial positions. We observe in figure [5|a) that 
the dipolar field (-82.2) is about twice as large as the component (i?2,6i) of field expulsion. In 
figure [5](b), the induction is maximum in the mid-plane z = 0, save for the (-B2,r) component 
for which the top and bottom boundary conditions play a major role. 



4. Higher orders 

As detailed in [16] , the interest of the iterative procedure is to associate induction effects 
with specific actions of velocity gradients. It is particularly convenient when a patterns 
develops through the iterations. As it is the case here, we define a tool to help us quantify 
the convergence of the pattern. Let Bj and B^ be the fields induced at respective orders j 
and k. We compute the scalar product 

(B.lBfc) = \— [ d^r B, ■ Bfc , (15) 

and the associated norm A/'(Bfc) = a/ (Bfe|Bfe). Comparisons are made using the normalized 
scalar product 

^ ^ , (16) 



Ar(B,) 



Bfc 
Ar(Bfc) 
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In the case of an azimuthal applied field, the expulsion eventually dominates. To wit, we 
compare in figure |6]the induced fields at order 1 and 3, and successive iterations of the axial 
field (B2fc). One finds Pi 3 ~ —0.9, while P2,3 ~ 0: the field induced at third order is almost 
exactly opposed to Bi and perpendicular to B2. As a result, the successive induction steps 
lead to the expulsion of the applied field, with Pk,k+i ~ and Pk,k+2 ~ —1, at higher orders. 
Bear in mind that this concerns normalized values. In dimensional units, one has : 

Bfc+2 ~ 7Bfc ,with 7 ~ -1/400 (17) 
After the lO*'^ order, we could not detect any appreciable evolution of the pattern. 



5. Evolution with Rm 

The magnetic Reynolds number is re-introduced in the summation Bind = YlT=i -^m-^fc- 
We have used the summed terms up to order 22, resulting in a 1% accuracy. We show in 
figure [7] the evolution with Rm of two components of the induction: the axial field in the 
center of the cylinder, and the toroidal field in the center of the columns. For Rm < 8, 
the calculation at second order yields a very good approximation of the net induction. 
This is interesting because the 2°'^ order truncation corresponds to the computation of the 
mean- field theory with a first order smoothing approximation [2D]- For magnetic Reynolds 
numbers greater than about 8, the contributions of higher orders in the summation need to 
be taken into account. One finds that the mean- field approximation tends to over-estimate 
the induced dipole field - figure [7](a), as well as the expulsion of the toroidal field - figure [7](b). 

Empirically, we observe that the radius of convergence of the integ er series is RJ^ — 17. 
This value can be understood from our observation that Pk,k+i ~ and Bk+2 ~ tB^. Indeed 
one can then rewrite the summation as 

CO 

B = {R^B, + RlB,) J2{-\^\Rl^f , (18) 

from which one immediately gets 1 / ^/[y\ ~ 20 for the radius of convergence. In addition, 
for Rm < R^, one gets: 

which shows that the divergence of the integer series actually lies in the existence of imagi- 
nary roots. Such a configurations is particularly suited to the use of Fade approximants [18] . 
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The result, plotted as a solid line in figure [7} shows that at large magnetic Reynolds numbers 
(R-m > 40) the axial induction and the expulsion may saturate. 



B. Induction from a radial field applied 

As detailed above, starting from an applied toroidal field, the a-effect generates a poloidal 
induced field with a large axial component, but it is easy to show that this axial component 
gives in turn very weak contributions to the induction. Indeed, an applied field in the axial 
direction induces fields that are two orders of magnitude weaker than the values obtained 
with other orientations. In contrast, we show in this section that a radially applied field 
generates an induced field which has a significant component in the azimuthal direction. 

Specifically, we consider an applied field Bq of the form 



Two remarks about this functional form are in order: 

(i) within the domain of resolution of the induction equation it is essential that Bq be 
divergence free because we solve iteratively AB^+i = —V x (v x B^) rather than AB^+i = 
-V VBfc - Bk ■ Vv. 

(ii) here, Bq is not divergence free for r < a, but in this domain the source term in 
AB 1 = — V X (v X Bq) vanishes with v. 

1. Induction at first order 

The induced field at first order is shown in figure [s] We observe in (a) and (b) that 
the magnetic induction is in phase with the flow in the columns, as opposed to the case 
of a toroidal applied field, where there was a phase difference of a quarter of a period. 
In the transverse plane, the effect of the azimuthal part of the velocity field is to twist 
the applied field lines, thereby generating azimuthal components in the magnetic field, in 
direct correspondence with the case of an azimuthal applied field, shown in figure |2](d). The 
magnetic field induced in the axial direction is due to the radial gradient of the vertical 
velocity and solves ABi z = —BQ{a/r)drV^. It corresponds to the stretching of applied field 
lines in the vertical direction. However, it is can be seen in figure |8ta,c) that the induction 




r 



(r < a = 0.1) 
(a < r) 



(20) 
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lies essentially near the boundary with the inner core (r ~ 0.6) compared to the boundary 
with the outer medium (r ~ 1). Since at both location there exists a large radial gradient 
in the vertical velocity, the difference has to be linked with the electromagnetic boundary 
conditions. The inner core has the same conductivity as the ffuid while the outer medium is 
insulating, so that the magnetic field is potential outside the cylinder. For a z-independent 
flow, one can show that this condition ensures that the induced field Bz vanishes at the 
r = 1 boundary [26j. In the case of our flow, the velocity field is not 2;-independent but its 
components have only second order variation with z in the mid-plane. One thus understands 
why the induction in the axial direction is so weak at r = 1, as shown in figures |8](a,c). 

2. Induction at second order 

As can be expected, in correspondence to the induction in the case of an azimuthal applied 
field, the screw motion in the columns produces again a current parallel to the applied field. 
This is shown in figure [oj^b) where one finds an induced radial e.m.f. maximum in the center 
of each column. This electromotive force has the symmetry of the helicity, reversing if either 
or are reversed (but not both simultaneously). The corresponding induced currents 
generate the axial field -82,2 shown in figure |9]^c). The field is periodic in the azimuthal 
direction with a period equal to half the period of the column distribution - in agreement 
with the Parker induction mechanism sketched in figure [sje (vertical arrows). 

When one computes the azimuthal average, a first observation is that the currents yield 
a non zero contribution (J2) which is essentially radial in the mid-plane z = 0. In addition 
the current lines close inside the conducting fluid, leading to the poloidal loops shown in 
figure |9](e) . They generate a field (B2) whose toroidal part is anti-symmetric about the z = 
plane, as shown in figure |9](f). 

One thus finds that a toroidal magnetic field is generated from the application of a radial 
magnetic field. The conversion is less efficient than the reverse process, discussed in the 
previous section. One computes max{(i?2,6i)/-Bo,r} ~ 7 ■ 10~^ instead of max{(_B2,r)/-So,0} ~ 
11 ■ 10^^. However, the major benefit is that the expulsion is much weaker. The component 
of (B2) opposed to the applied radial field is 10 times weaker than the induced toroidal 
field. This can be explained by the fact that, while B2 has a component opposed to the 
applied field in the center of the columns, it has a contribution which reinforces the applied 
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field in their periphery as can be seen in figure |9](d). Therefore, when the azimuthal average 
is taken, the expulsion effect is weakened. On the contrary, in the case of the azimuthal 
applied field, the expulsion effect from each columns add up collectively (see figure |3](b)), 
resulting in a larger contribution. 

Hence, expulsion is less effective in the direction perpendicular to the direction where 
scale separation develops (in the azimuthal direction the columns cross-section is an order of 
magnitude smaller than the cylinder diameter, while in the radial direction the characteristic 
size of the flow is equal to the width of the column). We conclude that in the case consid- 
ered here, scale separation does not particularly favor magnetic induction, but dramatically 
reduces expulsion. 



3. Higher orders 

The structure of the induced magnetic fleld is rapidly stabilized as higher orders are 
computed. As shown in flgure [TO] the flelds at orders 3 and 5 closely resembles that at order 
1; the normalized scalar product are Pi^s ~ —0.5 and Ps^s ~ —0.9. 

The flelds produced at a given order are again fairly orthogonal to the flelds at next or 
previous order : (B,fc|Bjt+i) ~ 0. In addition, as can be seen in flgure [Toj^g) , the iteration 
converges for even orders towards a quadrupolar structure with a negative feed-back in a 
two step mechanism: 

= -7Bfc (7 ~ 1/415) . (21) 



Evolution with Rn 



The Rm dependence of the induced radial and azimuthal magnetic fleld is shown in 



flgure 11, after summation of the terms up to order 22. Note that Be is sampled at (r = 
0.7, z = —0.6), because for a quadrupole, the toroidal fleld is very small in the z = plane. 
One observes in flgure [lT|(b) that a second order calculation is a correct approximation for 
magnetic Reynolds numbers up to 10. For higher Rm other terms need to be included. As 
in the case of the toroidal applied fleld, they tend to slow the increase of the induced fleld, 
mainly because of expulsion generated from rotational motion in the columns. The integer 
series diverges for Rj^ = 18, as expected from the value 1/ a/It] ~ 20. Results obtained using 
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the Pade approximants, although extending the computed induction beyond the radius of 
convergence of the series, do not point to a saturation at large Rm- 



IV. DYNAMO ACTION 



In the previous section, we showed that the induction mechanisms in the case of the 
T| flow consist in a mutual conversion between azimuthal and radial fields, through the 
a effect, along with an effect of expulsion of these fields by the rotating columns. We 
have identified two modes, dipole and quadrupole, mainly axisymmetric, that realize a 



feed-back loop in a two-step mechanism (equations (17) and (21)) with a negative sign, 
therefore leading to an anti-dynamo configuration. Following the ideas developed in [16], 
we can express these results using an induction operator formalism: from equation (|8]), we 
define £(Rm) = — Rm^"^ {V x (v x •)} for a velocity field v corresponding to the value 



Rm = 1. Equations (17) and (21) can then be interpreted in the following way: the (B2n) 



modes {n > 5) obtained in the induction studies in section III A (dipole mode) and IIIB 



(quadrupole mode) are eigenvectors of £^(Rm). From this observation, we will show in this 
section that using a poloidal/toroidal decomposition, a matrix analysis can be performed 
on the identified eigenmodes to find eigenvectors of £^(Rm) with positive eigenvalues, thus 
leading to possible dynamo solutions. 

Indeed, writing for simplicity C = £(1) (we then have £(Rm)B = Rm£B), let's assume 
that we can find a magnetic field Be that is not an eigenvector of the operator C, but 
of with a positive eigenvalue 7e. We can easily show that B^ = Bg + -^C^e is an 
eigenvector of C with a positive eigenvalue 7^ = Taking Rm = 1/^/%, we then have 

£(Rm)B5 = Bs, which defines B^ as a self-sustained magnetic field by the velocity field, at 
threshold Rm = 1/ 



A. Dynamo action in the flow. 

A dynamo cycle is often seen as a toroidal/poloidal feed-back loop [TT]. In the cylindrical 
geometry of the T| flow, it is easy to decompose an axisjTumetric field into its poloidal and 
toroidal part. We perform such a decomposition successively on the dipole and quadrupole 
axisymmetric eigenmodes of the operator and then compute the action of this operator on 
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the vector space generated by these two components. This leads us to a matrix formulation 
of the problem, allowing the identification of dynamo modes. We then perform the same 
analysis on a non axisymmetric mode. 



1. Generation of a dipole field. 



We start from the field (Bio) ^ figure |6](e) - which was shown (equation (17)) to be an 
eigenvector of C^, with a dipolar geometry and a negative eigenvalue. Figure 12 a) and (c) 
show the decomposition of this field into its toroidal - B J - and poloidal - B^ - components. 

Figure 12 b) and (d) show the axisymmetric part of the fields -C^Bj and C^B^- -^^ 



expected, the resulting fields are of opposite sign with respect to B^ and B^. In addition. 



their topology are very similar and look like a linear combination of the initial fields Bj and 



B^. Using the scalar products defined in equation (|15[), we can project (£^Bj) and (£^B^) 



onto the initial fields, thus defining an induction matrix: 

M.,r|) ^ f^'- '"^A . I ™ I ) 1^-1 , (22) 

\Mtp Mtt) {({C'B^) I Bf) ({C^BJ) I BT 

where we have taken A/'(Bj) = A/'(B^) = 1, using the norm Af defined in section III A 4 
This matrix Mrf(Ti) is the restriction of the two-step induction operator to the vector 

space of the axisymmetric dipoles. The diagonal terms represent the expulsion effect and 

the extra-diagonal terms represent the action of the a effect. 

In the case of the flow with {d/r = 0.4, ^ = 1.25 and = 4), we compute: 



Mrf(Ti) = -10-^ (23) 




This matrix has a negative eigenvalue = —25.10"'^ = —1/400, corresponding to the 
values obtained by iterating the induction operator in the case of an azimuthal or radial 
applied field. In addition, it has a positive eigenvalue = 10~^, showing that the flow can 
actually sustain an axial dipole dynamo, for > = 1/a/aJ = 100. The eigenmodes 
corresponding to both eigenvalues have the following structure: 
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B+ = 0.88B^ - 0.47BJ 



B, =0.81B^ + 0.58BJ 



(24) 



We have assumed so far that (£^Bj) and (£^B^) belong to the vector space generated by 
BJ and B^. In order to check the vahdity of this hypothesis, we show some comparisons in 



figure 13, using profiles normalized by their maximum value: in (a) the axial profiles of the 
radial component of (£^Bj) and (£^B^) are compared to the equivalent profile for B^. One 
can see that the overlap is very good. In (b), the radial profiles of the azimuthal component 
of (£^BJ) and (£^B^) are compared to the radial profile for Bj. In the case of the a-effect 
((£^B^) profile), the overlap is very good again. In the case of the expulsion effect ((£^Bj) 
profile), there is a discrepancy. We have computed the error to be about 10%. Thus the 
Mtt term in the matrix represents the expulsion mechanism for the applied toroidal field 
within an error of 10 %, while the other elements of the matrix can be considered as correct 
within less than 1 %. 



2. Generation of a quadrupole field. 

We follow the same procedure as for the dipole case. We start from (-Bio) presented in 
figure QO[g), which was shown (equation (21)) to be an eigenvector of the operator C^, with 



a quadrupolar geometry and a negative eigenvalue. This mode is orthogonal to the dipole 
mode used in the previous section. We split this field into its toroidal - B^ - and poloidal - 



Bi^ - components (figure 14 a) and (c)). Figures 14 b) and (d) show the axisymmetric part 



of the fields >C^B^ and >C^B^. Again, the resulting fields are of opposite sign with respect to 
B^ and B^ and have very similar topology. We define the matrix M as in the dipole case, 
yielding: 



M,(ri) = -10- 



10 12 



12 11, 

This matrix has two real eigenvalues of opposite sign: = 2.10"^ and 
-1/430. The flow can thus also sustain an axial quadrupole dynamo, with a threshold 



(25) 

-23.10^^ = 
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= 1/ = 80. The two corresponding eigenmodes have the following structure: 

B+ = 0.72B^ - 0.69B^ 

(26) 

Eg = 0.69B^ + 0.72BJ 
As in the dipole case, the error made assuming that (B2) belongs to the vector space 
generated by B^ and B^ can be estimated by comparing the profiles of (£^Bf) and {C'^'B'^) 
to B^ and B^. We find here also that the only component that does not overlap correctly 
is the one corresponding to the expulsion of the azimuthal applied field - M^t ^ with an 
error of about 10%. 



3. Transverse dipole. 

Until now, we have only considered axisymmetric fields. However, since the fiow 
presents several analogies with the Roberts fiow, it would be interesting to study the possi- 
bility of generating a transverse dipole (perpendicular to the columns), as observed in the 
Karlsruhe dynamo [H]. We follow the same strategy as before: a uniform transverse field 
Bq = B^e^ is applied to the fiow and we compute the fields B^ obtained after k iterations. 
They rapidly converge towards a stable structure. After 8 iterations, the fields B^ and C^B^ 



have an overlap close to 100%, and the eigenvalue is negative, 7 = —1/400. Figure 15 shows 
two cuts of the field Bg. 

In order to use the same matrix analysis as in the axisymmetric cases, one must find a 
way to decompose Bg into a couple of vectors generating a subspace closed under No 
poloidal/toroidal decomposition can be made here since the mode is not axisymmetric, but 
symmetry considerations can be used, based on the following observation: in section IV A 1 



and IV A 2 the toroidal/poloidal decomposition of the eigenvector also resulted in separating 
the symmetric and antisymmetric parts of the eigenvector under the refiection symmetry 
with respect to the plane z = 0. In what follows, we will call Sz this symmetry. One can 



then notice in figure 15 that the x component of the field is symmetric under S^, while the y 
component is antisymmetric. The system being finite and the exterior being insulating, the 
X and y components are each associated to an axial component having the same symmetry: 
Bx, symmetric with respect to Sz, is associated with a symmetric axial component, i?f ; By, 
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antisymmetric, is associated with the antisymmetric part of the axial component, -B". We 
thus define the vectors B° = -B^e^. + -B^e^ and B° = ByBy + B'^e^ as the generating vectors. 

Applying to these vectors, one observes that, in this case too, the resulting vectors 
can be expressed as linear combinations of and B°. can therefore be expressed in this 
basis as the matrix : 

M,(Ti) = -10-^ I (27) 

This matrix has a positive eigenvalue = 1.6 10~^, showing that the T| fiow is also 
able to sustain a transverse dipole, with a threshold equal to 79, very close to the one of the 
axial quadrupole. 

In a fiow, we thus observed that several dynamo modes can be sustained when Rm — 
80. This result is consistent with other numerical studies in thermal convection P2], which 
have shown that various dynamo solutions can coexist in the same region of the parameter 
space. This kind of behavior has also been observed in the VKS dynamo experiment [35]. 



B. Dynamo mechanisms in the Tl flow. 



1. Symmetry considerations. 

The analytical expression of the Tl fiow is given in equation Q. Compared to the T| 
fiow (equation ([s])), the only difference is that the azimuthal and axial components are now 
respectively symmetric and antisymmetric under Sz- Tl is a superposition of two Tl fiows, 
one in each half-cylinder, symmetric with respect to Sz- They have opposite helicity, since 
their axial component is reversed, while the rotation of the columns is unchanged. The 
strong similarity between the T^ and Tl fiows indicates that the same mechanisms, a effect 



and expulsion, will take place. Figure 16 compares the schematic induction mechanisms for 
the a effect in both types of flow. In (a), the case of the T^ flow is recalled, where the 
induced current (j = cr«B) is parallel to the applied fleld with the same sign, since helicity 
Ti. is negative in this flow and the a coefficient is proportional to — 7i. In (b), it is shown 
that in the Tl flow, applying a symmetric toroidal fleld results in a symmetric poloidal 
field, corresponding to a quadrupolar geometry. In (c), we show that the dipolar geometry 
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(antisymmetric poloidal field) is obtained by applying an antisymmetric toroidal field. 



2. Search of axisymmetric dynamo modes. 
Dipoleir mode: 

The symmetry considerations of the previous paragraph suggest that a dipolar mode 
can be obtained in the Tl fiow by an iterative application of C to an initial toroidal field 
antisymmetric with respect to Sz- We chose an initial field of the form Bq = sin(7r^)e6i. 
The iterations converge rapidly towards a mode whose toroidal and poloidal parts, once 
the azimuthal average is done, form a basis which is closed under the action of The 
corresponding matrix is: 



It has two negative eigenvalues Xf — —24 10~^ and A2 = —1.2 10~^, so that with the 
current parameters {d/r — 0.4, ^ = 1.25 and n — A), the Tl fiow is unable to sustain an 
axial dipole dynamo mode. 

Quadrupolair mode: 

In this case, the symmetry study shows that the initial azimuthal field must be symmetric 
under Sz- Two choices can be tested: an azimuthal ^-independent field or a radial field 
of the form Bq = Bq^ cos{TTz/H)er. In both cases the iterations converge towards fields 
with identical structures. As previously, the toroidal and poloidal components of the field 
obtained after 10 iterations define the basis for the quadrupole vector space. In this basis, 
the expression of C"^ is: 




(28) 




(29) 



Again this matrix has two negative eigenvalues 



-20.10-4 and A| = -0.2.10"^ 
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It is interesting to note that in all the cases where a dynamo was possible, the 
corresponding matrix (Mrf(Ti), Mg(T^) and Mj(Ti)) presented the characteristics that 
the product of the non-diagonal terms was larger than the product of the diagonal 
terms. One can easily show that this is a necessary and sufficient condition for a 2x2 
matrix with negative diagonal coefficients in order to have a positive eigenvalue. And 
this condition physically corresponds to the facts that the expulsion mechanism (diagonal 
terms) is weaker than the a mechanism (non-diagonal terms). On the contrary, in the 
case of (Mrf(Tl)) and (Mg(Tl)), this condition is not met. In the quadrupolar case, both 
products are very close, resulting in a small value for \Xl\ (close to 10~^). We can therefore 
expect that by increasing the C, parameter (reducing the rotational component of the flow, 
responsible for the expulsion effect), a quadrupolar dynamo could be observed in the Tl flow. 



Transverse mode: 



In the case of the flow, we have iterated the C operator starting from a uniform 
transverse field and then defined the basis vectors using as a criterion the behavior of the 
cartesian components of the converged mode under the Sz symmetry. We tried to follow 
the same procedure in the case of the Tl fiow, but it turned out that the converged mode 
is formed of two dipolar structures at 90 degrees from each other, having the same behavior 
with respect to the symmetry. Therefore, it is not possible to use this symmetry to 
construct the basis vectors. And indeed, we were not able to find any basis that would be 
closed under the action of C^. This shows a limit of our method. 



C. Generalization for other values of the poloidal/toroidal ratio. 

1. General expression for the matrix M. 

Until now, all our studies have been based on the same fiow parameters [d/r = 0.4, 



^ = 1.25 and n = 4). As we noticed in section IV B 2 that a change in the axial/rotational 
ratio parameter) might allow the Tl fiow to sustain a quadrupolar dynamo, we redefine 
the C operator as a linear function of the poloidal and toroidal components of the fiow. We 
normalize the components v"^ and of the velocity field: max(v"^) = 1 and max(v'^) = 1 
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and rewrite equation as V = Va'v^ + Vr\^ where Va and Vr represent the maximum 
amphtude of each component. Let Ca be the induction operator when V = v"^ and Cr be 
the induction operator when V = v^; C can then be written as the hnear combination: 



/: = VACA + VnCR , (30) 



yielding: 



= V^CajCa + VaVr{CrCa + CaCr) + V^CrCr , (31) 
so that the 4 matrix elements Mjj can also be written as quadratic forms of Va and Vr. 

= aijVl + hjVRVA + Ci.Vl (32) 

As we studied the induction mechanisms, we noticed that the a effect and the expulsion 
behave differently under a reversal of the axial pumping {Va —>■ —Va) or of the columns 
rotation (Vr —Vr). More precisely, it was observed that the expulsion effect is indepen- 
dent of these sign changes, whereas the a effect transforms as the product VaVr. These 



observations allow to eliminate some terms in equation (32), yielding: 



I aV^ + bVl cVrVa , 
M{Vr,Va)=[ ^ ^ 33 
dVRVA eV^ + fVl' 

2. Positive eigenvalue diagram for the flow. 

In the case of the T| flow, the value of these 6 coefficients are computed using successively 
the dipole basis, the quadrupole basis and the transverse basis obtained respectively in 



sections IV A 1 IV A 2 and IV A 3 The corresponding expressions for the matrix M are 



the following: 
Axial dipole: 



, 22Vi + 6VJ 46VaVr 



M,iVA,VR) = -10-' ^ . , \J (34) 

ISVaVr 14(V| + VI) 
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Axial quadrupole: 



, , 23V? + 2WaVr , 

M,(V4,V^) = -10-M ^ ^ _ . I (35) 



Transverse dipole: 



2WaVr lliYl + V^) 



_4 , 22V:| + 13Vj 291/^ Vr 



Mt(VA,VR) = -10-^ " " (36) 

46v:iyR i6v;| + uv^^ 

From these expressions, the largest eigenvalue Xmax of each matrix can be computed as a 



function of Va and Vr. The result is shown in figure 17 In order to facilitate the reading of 
these plots, when Xmax was negative (no dynamo), we artificially set its value to zero. For 
each case, two regions are evidenced: the first one, for which Xmax > 0, corresponds to the 
possibility to observe a dynamo for the considered geometric parameters { d/r = 0.4 and 
n = 4) with a threshold = l/V-^max- The second, for which Xmax < 0, corresponds to 
the case where no dynamo instability can take place. Figures (a), (b) and (c) show that the 
three dynamo modes coexist in the same region, in the neighborhood of the line Va = Vr. 
This is not surprising, since when one of the velocity component dominates the other, the 
expulsion mechanism is more important than the a effect, which needs both components 
together. 

For a given couple (Va,Vr), one can observe that Xq is always larger than (corre- 
sponding to a lower threshold). In the same way, the transverse dipole always has a lower 
threshold than the axial dipole. A summary of the predominance of the different modes is 
given in figure [T7|^d), where one can see that the transverse dipole is favored by higher ratio 
axial/rotational, whereas the quadrupole is favored by a lower ratio. 

3. Positive eigenvalue diagram for the Tl flow. 

We performed the same analysis for the Tl fiow. The matrix M has the following ex- 
pressions: 

Axial dipole: 

^ ,' 16V1 + 8V^ 33VaVr , 
M^iVA, Vr) = -10~' \ ^ ^ ^ ^ (37) 

IGVaVr 12{VI + VJ) 
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Axial quadrupole: 



M^{VA.Vn) = \ ^ ^ , J (38) 

2WaVr llVl + 8V| 



Figure 18 shows tlie evolution of \max in the case of the quadrupole mode. It confirms the 



assumption made in section 1VB2, namely that starting from the value ^ = 1.25 studied 



in that section (cross in figure 18), one could evolve the system towards a dynamo state 

to a 



by reducing the intensity of the rotation (increasing ^, corresponding in figure 18 
displacement towards the top left corner). 

On the other hand, there is no value of the couple (V^, Vr) for which the axial dipole can 
be sustained by the Tl flow. Some studies have shown [HTj that in order for a Tl-type flow 
to sustain an axial dipole, the presence of differential rotation, absent in our study, would 
be necessary. 



V. CONCLUDING REMARKS. 



A better understanding of the MHD induction mechanisms in a given system can 
help to build dynamo cycles. In the first part of this study, we have identified two 
mechanisms, related respectively to the a-effect and to expulsion by vortices. We expressed 
these mechanisms in terms of the induction operator and, using a poloidal/toroidal 
decomposition of the eigenvectors of we were able to perform a matrix analysis 
leading to the determination of self-sustained magnetic modes. The benefit is to give a 
complete description of the modes whose interaction lead to dynamo action. In regards 
to natural or experimental conditions, it helps understand which features of the velocity 
field favor or hinder dynamo action. In addition, as already noticed in other studies [28] , 
the competition between the a-effect, favorable to the dynamo, and the expulsion effect, 
that works against it, can be monitored by the ratio of poloidal to toroidal components 
of the velocity field (in our case the axial to rotational ratio Va/Vr). Our studies show 
that a positive feed-back requires a comparable amplitude for rotation and pumping. In 
the case of the flow (for 8 columns and an aspect ratio of 0.4), we observed that both 
axisymmetric modes (dipolar and quadrupolar) can be sustained, with a threshold of the 
order of Rm = 100. We also showed that a transverse dipole mode can exist, as one 
could expect, because of the analogy between our flow and the Karlsruhe dynamo. On 
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the other hand, the Tl flow can only sustain a quadrupolar mode, when the pumping 
amphtude is larger than the rotation amplitude. For the axial dipole mode, the poloidal to 
toroidal conversion seems too weak to compensate the strong expulsion of the azimuthal field. 

Coming back to the Earth's case, our model system leads to several observations 
which may be relevant. An dynamo process relies on the helicity contained in Busse's 
columns but Eckman pumping would give a very weak source of axial motion, since the 
Eckman number E is of the order of 10~^^ and the ratio of the axial flow to the rotational 
flow scales like E^^'^ [30] • Another source of axial velocity could be by the /9-effect due 
to the curvature of the core-mantle boundary - note that in this case the axial flow is in 
phase with the radial flow rather than with the vorticity [31]. A large scale dipole fleld 
could also be generated from an a — dynamo. It would require differential rotation 
as provided, for instance, by zonal winds [331 El] or super-rotation effects as observed 
in a recent laboratory experiment [32] These ingredients could in principle be added 
to the model studied here, and the procedure used to determine which dynamo modes 
(dipole, quadrupole or other) are likely to exist for a given range of Reynolds numbers. As 
proposed for geomagnetism [36], and recently observed in the VKS experiment [35j, the 
close proximity of dynamo modes may be essential for the development of dynamical regimes. 
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Volk et al, Figure 1 




FIG. 1: Geometry investigated, (a) Sketch of the column arrangement for the flow, (b) 
Column arrangement for the Tl flow, (b) Cut of the flow in the plane z = (colors correspond 
to the vertical flow, arrows correspond to the flow in the plane z = 0) (d) Cut of the flow in 
the plane ^ = 0. (e) Cut of the Tl flow in the plane 9 = 0. 
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Volk et al, Figure 2 
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FIG. 2: T^: magnetic field induced at first order for an applied toroidal field -Bo.e- Fields cuts 
are shown in the z = mid-plane, (a) induced field along the axial direction, Bz', (b) transverse 
components {Br,Bg). Associated distortion of the applied field lines: (c) axial induced field from 
axial stretching between columns, (d) radial induced field from the rotational flow. 
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Volk et al, Figure 3 





FIG. 3: T^: magnetic field induced at second order for an applied toroidal field Bq q. (a): vertical 
component of the electromotive force 62,2, shown in the mid plane z = 0, (b): transverse compo- 
nents (e2,r;e2,6») in the same z = plane; (c) corresponding induced magnetic field i?2,2 and (d): 
(i?2,rj B2fi)\ (e): schematics of the cooperative induction effect from a pair of neighboring columns, 
generation of -62,2; (f): associated expulsion mechanism, generation of B2fi. 
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Volk et al, Figure 4 
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FIG. 4: T^: azimuthal average of the fields induced at second order for an applied toroidal field 
Bo g. Azimuthal averages of the fields are shown in an (r, z) plane, with poloidal components 
indicated by arrows and azimuthal component in color scale, (a): current (J2); (b): corresponding 
magnetic field (B2). 
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Volk et al, Figure 5 
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FIG. 5: T|: spatial profiles of (B2). (a): variation with r in the z = plane; (b): axial variation 
at r = 0.8 (the center of the columns). The magnitude of the applied field is 0.25 in the center of 
the columns. 
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Volk et al, Figure 6 
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FIG. 6: magnetic field induced at higher orders for an applied toroidal field Bq q. (a,b): axial 
and transverse components of Bi; (c,d): corresponding plots for B3; (e): evolution of the azimuthal 
average ((B2), (B2fc), . . . (Bio)) 
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Volk et al, Figure 7 
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FIG. 7: T^: evolution with Rm. (a): axial field in the center of the annulus, at {r = 0, z = 0); 
(b): toroidal field in the columns, at (r = 0.8,2: = 0). RJ^ is the radius of convergence of the 
integer series. The dashed line corresponds to the summation stopped at order 2, the continuous 
line with dots to the summation up to order 22 and the continuous line to the summation using 
Pade approximants. 
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Volk et al, Figure 8 
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FIG. 8: T^: magnetic field induced at first order for an applied radial field -Bo,r- Field cuts are 
shown in the z = mid-plane, (a) induced field along the axial direction, Bz', (b) transverse 
components {Br,BQ). (c) Radial profiles of the induced fields, at (z = 0,0 = 0). 
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Volk et al, Figure 9 
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FIG. 9: T|: induction at second order for an applied radial field -Bo,t-- (a-jb) electromotive force; 
(c,d) induced field B2; (e,f) meridian view, for the azimuthally averaged fields (J2) and (B2). 
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FIG. 10: T^: induction at higher orders for an apphed radial field -Bo,r- (a-f): structure of the 
field induced at orders 1, 3 and 5, shown in the mid plane z = 0. (g): structure of the azimuthal 
averaged (J„) of the fields induced at orders n = 2,4, 6, 8, 10. 
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FIG. 11: T^: induction at higher orders for an apphed radial field i?o,r- Variation with Rm of the 
induced fields, (a) radial induced component, sampled at (r = 0.7, z = 0). (b) azimuthal induced 
component, sampled at (r = 0.7, z = —0.6). The dashed line corresponds to the summation stopped 
at order 2, the continuous line with dots to the summation up to order 22 and the continuous line 
to the summation using Pade approximants. 
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FIG. 12: Dipolar mode: (a),(c) resp. toroidal and poloidal component of the dipole eigenvector, 
(b), (d), axisymmetric part of the fields obtained by applying C? to the vector fields (a) and (c). 
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FIG. 13: Dipolar mode: estimate of the error made doing the projection (a): axial profile of Bj. at 
r = 0.8 from the Oz axis, (b): radial profile of Bg in the median plane. (■) : Components of the 

») Components of (£^BJ). 



applied field (♦) Components of (£^B^) 
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(C) Bf (d) {C^bP) 

FIG. 14: Quadrupolar mode: (a),(c) resp. toroidal and poloidal component of the quadrupole 
eigenvector, (b), (d), axisymmetric part of the fields obtained by applying to the vector fields 
(a) and (c). 
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FIG. 15: Transverse dipole mode: (a) {Ox,Oz) cut. (b) {Oy,Oz) cut. 
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FIG. 16: Comparison between the schematic induction mechanisms for the and Tl flows, due 
to the a efi'ect in the case of a toroidal applied field. The black line represents the applied field, 
the green dashed line the resulting order 2 current and the red line the resulting order 2 magnetic 
field, (a) fiow (b) Tl fiow, with a symmetric applied field (c) Tl fiow, with an antisymmetric 
applied field. 
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FIG. 17: flow : evolution of the largest eigenvalue of the induction matrix in the {Vr,Va) 
plane, (a) axial dipole. (b) axial quadrupole (c) transverse dipole (d) predominance diagram for 
the different possible modes. For simplicity, when the eigenvalue is negative we plotted a zero 
value. On each plot, the cross marks the parameter couple (Vr,V^)=(0.83,0.53) of the particular 
flow studied in IIV A[ 
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FIG. 18: Tl flow : evolution of the largest eigenvalue of the induction matrix in the {Vr,Va) 
plane, for the axial quadrupole. For simplicity, when the eigenvalue is negative we plotted a zero 
value. The cross marks the parameter couple (Vr,,Va)=(0.83,0.53) of the particular flow studied in 

HYll 
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